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Abstract 

The SafFman- Taylor problem addresses the morphological instability of an 
interface separating two immiscible, viscous fluids when they move in a nar- 
row gap between two flat parallel plates (Hele-Shaw cell). In this work, we 
extend the classic SafFman- Taylor situation, by considering the flow between 
two curved, closely spaced, concentric spheres (spherical Hele-Shaw cell). We 
derive the mode-coupling differential equation for the interface perturbation 
amplitudes and study both linear and nonlinear flow regimes. The effect of 
the spherical cell (positive) spatial curvature on the shape of the interfacial 
patterns is investigated. We show that stability properties of the fluid-fluid 
interface are sensitive to the curvature of the surface. In particular, it is found 
that positive spatial curvature inhibits finger tip-splitting. Hele-Shaw flow on 
weakly negative, curved surfaces is briefly discussed. 



Typeset using REVTgX 



1 



I. INTRODUCTION 



Formation and evolution of dynamic structures is an exciting area of nonlinear phe- 
nomenology. Of particular practical and theoretical interest is the hydrodynamic pattern 
formation of the growing interface between two fluids. One of the best studied pattern- 
formation systems of this type is the Saff man- Taylor problem it addresses motion of 
two viscous, immiscible fluids in the narrow space between two parallel, flat plates known 
as a Hele-Shaw cell. When a fluid of low viscosity displaces a fluid of higher viscosity, the 
interface between them becomes unstable and deforms. Dynamic competition leads to the 
formation of fingering structures. 

Experiments and theory focus on two principal geometries: (i) rectangular [|l],H and (ii) 
radial In rectangular cells the unperturbed interface is straight and the unperturbed 

flow is uniform and parallel to cell walls. In the radial case the unperturbed interface is 
circular with the less viscous fluid pumped into the more viscous one at a point and the 
flow radially outward. For both situations, the initial development of the interface instability 
tracks the predictions of linear stability theory . After the initial surface deformation, as 
the unstable modes of perturbation grow, they become coupled in a weakly nonlinear stage 
of evolution Finally, the system evolves to a complicated late stage, characterized 

by formation of fingering structures, in which nonlinear effects dominate p|. As a result, 
beautiful finger-bubble undulated structures are formed in rectangular cells, while visually 
striking, fan-like, branched patterns rise in radial flow. Spreading, splitting and competition 
are the three basic growth mechanisms of the viscous fingering process 0|. In this work we 
are particularly interested in tip-splitting events. 

Despite the extensive experimental and theoretical work on the Saffman- Taylor problem 
in both radial and rectangular setups, the majority of the studies focus on flow in flat 
Hele-Shaw cells. Curiously, the dynamic behavior for flow in flat cells is described by the 
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very same set of equations as those for flow in porous media 0, which is indeed a highly 
non-flat environment, characterized by voids and curved internal surfaces. Even though 
viscous fingering is not restricted to occur between flat surfaces, the study of the Saffman- 
Taylor problem on curved surfaces has been practically neglected. One exception is the 
experimental work by Zhao and Maher which considers flow in a cylindrical Hele-Shaw 
cell with a large radius of curvature. The flow is performed parallel to the cell axis. They 
used the cylindrical cell as an experimental realization of a flat rectangular Hele-Shaw cell 
with periodic boundary conditions. Their experiments showed that instabilities in cells with 
periodic boundary conditions are qualitatively similar to instabilities in cells with physical 
sidewalls. The authors did not explore the influence of cell curvature. 

On the theoretical side, Entov and Etingof |^ considered the general mathematical prob- 
lem of viscous flow in non-planar Hele-Shaw cells, in the zero-surface-tension limit. They 
applied conformal mapping techniques to derive a class of explicit solutions for the shape of 
the fluid domain under study. As in reference 0, the authors in were not interested in 
the influence of the cell curvature in the shape of the interfacial patterns. 

Spatial curvature has been a relatively overlooked feature in the study of other pattern 
formation systems as well. Only recently researchers started to investigate how the properties 
of the patterns could be affected by the curvature of the surface in which such structures 
evolve [jlO|-0. In the early 1990's Levine and collaborators studied coarsening of two- 



dimensional foams [IC] and grain growth on curved surfaces. They showed that the 
stability properties of such froth bubbles and grains depended on the curvature of the surface. 



Ref. [|T2[ reviews the dependence of many chemical and biochemical surface processes on 
curvature. Numerical studies of reaction-diffusion systems in curved spaces examine the 



evolution of spiral waves [13| and the occurence of Turing patterns on a sphere [14]. These 



numerical simulations |13,14| indicate that curvature imposes geometrical restrictions on the 
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shape of the patterns. Finally, Schoenborn and Desai [|T5| studied the intra-surface kinetics 
of phase ordering on curved surfaces. 

For the Saffman- Taylor problem, the interplay between Hele-Shaw cell curvature and 
interfacial pattern formation is largely unexplored in the present literature. However, in a 
recent mode-coupling analysis of radial flow in flat cells, Miranda and Widom |^ suggested 
that cell curvature could be used as a control parameter to regulate the tendency towards 
finger tip-splitting. A thorough investigation of the relationship between cell curvature and 
the fiuid-fiuid interface dynamics still needs to be addressed. In this work we begin such 
investigations, focusing on the Saffman- Taylor problem on a sphere. Gravity effects have 



been studied in reference [0. 

The study of viscous flow in a nonplanar Hele-Shaw cell is of interest for both scientific 
and practical reasons. On the scientific level, the influence of spatial curvature on hydrody- 
namic flow is a matter of fundamental interest. It also provides a simple mathematical model 
to describe more general situations involving the filling of a thin cavity between two walls 
of a given shape with fluid. On the practical level, it may have applications in a number 
of industrial, manufacturing processes, ranging through pressure moulding of molten metals 
and polymer materials and formation of coating defects in drying paint thin films |]18 



The outline of the work is the following: section || defines Hele-Shaw flow between con- 
centric, thin spherical shells, considering fluid injection (withdrawal) at the north (south) 
pole. We derive a differential equation describing the early nonlinear evolution of the inter- 
face modes. In section |nT| , we interpret results obtained in section || and investigate both 
the linear and weakly nonlinear evolution of the system. The flat space limit of infinite 
radius of curvature is examined in section |111 A| . Section lllB discusses the linear growth 
rates of unstable modes and relates these to the degree of lattitude of the unperturbed 
interface. Nonlinear analysis in section pTO concentrates on the effect of cell curvature 
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on finger tip-splitting. We show that positive spatial curvature inhibits splitting. Flow on 
the northern and southern hemispheres are contrasted, and a symmetry-breaking behavior 
is detected: for the southern hemisphere tip-splitting is replaced by finger tip-sharpening. 
Flow on weakly negative curved surfaces is briefly discussed. Section |^ presents our final 
remarks. An appendix derives Darcy's law for flow between concentric spheres. 



II. THE MODE COUPLING DIFFERENTIAL EQUATION 

Consider two immiscible, incompressible, viscous fluids, flowing in a narrow gap of thick- 
ness b, between two concentric, thin spherical shells (see flgure 1). We name this device 
the spherical Hele-Shaw cell. Assume that b is smaller than any other length scale in the 
problem, so the system is effectively two-dimensional. The flow takes place in the surface 
of a two-dimensional sphere, embedded in three-dimensional space, and endowed with the 
metric 

ds' = dp' + a^sm^(^^^dif'', (1) 

where a is the radius of curvature of the sphere, < < 27r denotes the polar angle 
measured on the sphere and < p < vra is the geodesic distance from the radial flow source, 
located at the sphere's north pole. The sphere has a constant, positive Gaussian curvature 
K = l/a^ and a constant mean curvature H = 1/a. We presume the Gaussian curvature is 
more relevant than the mean curvature for reasons discussed in sections 1111 A| and [111 CI. At 



any stage of our calculation, the "flat-cell" limit a ^ oo (or, equivalently K ^ 0) gives all 
the well-known familiar results for flow in flat Hele-Shaw cells. 

Denote the viscosities of the upper and lower fluids, respectively as rji and ri2. Between 
the two fluids there exists a surface tension a (flgure 1). The flows are assumed to be 
irrotational, except at the interface. Fluid 1 is injected into fluid 2 through an inlet located 
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at the sphere's north pole, at a given flow rate Q, which is the area covered per unit time. 
Fluid 2 is simultaneously withdrawn, at the same rate, through an outlet placed at the south 
pole. 

During the flow, the fluid-fluid interface has a perturbed shape described as p = 71 = R+ 
({(p, t). The interface perturbation amplitude is represented by ({ip, t), and R = R{t) denotes 
the time-dependent unperturbed radius. We can calculate R{t) from the time- dependent 
surface area A(t) = Ana? sin^(i?(t)/2a). For flow rate Q we write A{t) = Ana? sin^(/2o/2a) + 
Qt then solve for 

R{t) = a arccos (^Cq - ^^^^ , (2) 

where Cq = cos{Ro/a), and Rq is the unperturbed radius at t = 0. The unperturbed shape 
is a polar cap of geodesic radius p = R, surface area A and circumference C = 27ra sin (i?/ a). 
Note the identity Q = vC where v = dR/dt is the velocity of the unperturbed interface. 
We express the net perturbation C{ip,t) in the form of a Fourier expansion 

C{v,t) = J2 Cn{t)exp{inip), (3) 

72= — CO 

where 

Cnit) = — / C(¥',^)exp(-m¥?) dip (4) 
Arc Jo 



denotes the complex Fourier mode amplitudes and n=0, ±1, ±2, ... is the discrete azimuthal 
wave number. To keep the area of the perturbed shape independent of the perturbation (, 
we express the Fourier mode corresponding to n = in the expansion as 

Co(t) = -^cot(^)^|C„(t)p. (5) 



The constraint is intrinsically a nonlinear concern and is not required in linear stability 
analysis. 



Since we are interested in both linear and early nonlinear behavior of the system, we 
must derive a differential equation for correct to second order. This second-order mode- 
coupling equation considers the presence of a full spectrum of modes. We begin the hy- 
drodynamic study of the system by considering a generalized version of the usual Darcy's 
law ll|J^, adjusted to describe flow between concentric spheres (see Appendix ^ 

6^ 
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Vp„ (6) 



where = \-j{p,(p) and pj = pj{p,ip) are, respectively, the velocity and pressure in fluids 
j = I and 2. The gradient in equation (H), associated with the metric (|l|), is [p^j 

^ d ^ 1 d ^ 

op asm[p/a) oip 

where the unit vectors p and point in the direction of increase of p and ip, respectively. 
Equation (|^ was obtained from the corresponding three-dimensional expression for the 
gradient in spherical coordinates (r, 6, by keeping r = a and noting that 6 = p/a. 

At the interface, the pressure difference between the two fluids is governed by the mean 
curvature (the average of the two principal curvatures) of the fluid- fluid interface 0. We 
can identify the directions of the principal curvatures in the limit of b smaller than any 
other length scale by considering a "tangent Hele-Shaw cell" consisting of two parallel planes 
tangent, respectively to the inner and outer spheres at some point along the interface. Within 
the tangent Hele-Shaw cell, one principal curvature is associated with the interface profile 
in the direction perpendicular to the tangent planes. We call this curvature k± and note 
that it is of order 1/6 with a specific value set by interface contact angles. The remaining 
direction of principal curvature is parallel to the tangent plane and tangent to the interface. 
We call this curvature k\\. Then the pressure jump boundary condition may be written 

{Pi - P2)\-R = {i^W + i^±) In- (8) 
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As was the case for flow in flat Hele-Shaw cells, k± is much larger than /ty but is is nearly 
constant [2^,21]. This curvature does not significantly affect the motion in our problem, 
because its gradient is nearly zero. 

Since the closed boundary describing the fluid-fluid interface is itself on the top of a 
curved surface (sphere), the calculation of the "intra-surface" interface curvature k\\ is not 
as simple as it was in the flat-cell case P|,p!9|. Taking into consideration the fact that the 
interface evolves in the surface of a sphere of radius a, we derive a slightly involved expression 
for the fluid-fluid interface curvature 

cos{p/a) a^sin^{p/a) +2cos{p/a) - a sin(p/a) 

i. (9) 



a2sin2(p/a) + (f 



3/2 



The sign convention for the curvature is such that a circular interface above the equator 
has positive curvature, whereas it has negative curvature below the equator. Keeping terms 
up to second order in the perturbation amplitude C, we rewrite the interfacial curvature as 



C 



aS a^S"^ 



1 / d^C\ c 



a?S^ 



.2 1 (dCV 



(10) 



where we have introduced the shorthand notation 5* = sm{R/a) and C = cos{R/a). Exactly 
at the equator (C = 0), only the term in ([l0|) that is linear in ( survives. 

Taking advantage of the irrotational and incompressible flow conditions, we define the 
velocity potential (pj in each of the fluids, where Vj = — V0j. The velocity potential satisfies 
Laplace's equation V^0j = 0, where is the two-dimensional Laplacian defined on the 
surface of the sphere. Combining the velocity potential with equations (||) and (H) for the 
pressure difference and the generalized Darcy's law (BI), we write the equation of motion 
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where 



A 



V2 - Vi 
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(12) 



is the viscosity contrast and 
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a. 



(13) 



12(771 + 772) ' 
contains the surface tension. 

Now define Fourier expansions for the velocity potentials 0^. Far from the interface the 
velocity field should approach the unperturbed steady fiow with a circular interface of radius 
R. Thus for p ^ and p — > vra the velocity potentials 0j approach 0°, the velocity potentials 
for purely radial (p direction) fiow, satisfying Laplace's equation 

tan(p/2a) 



Q , 

log 

27r ^ 



J' 



(14) 



_tan(i?/2a) 

where are independent of p and The general velocity potentials obeying all these 
requirements are 



tan(i?/2a) 



(-ly \n\ 



exp(m(y9). 



(15) 



tan(p/2a) 

The trigonometric dependence on p transforms, in the fiat-cell limit a ^ cxd, into the ratio 
Rj p. In order to calculate the mode coupling differential equation for the system, we sub- 
stitute expansions ([T^) and ([T5|) into the equation of motion (|TT]), keep second order terms 
in the perturbation amplitudes, and Fourier transform them. 

To conclude our derivation we need additional relations expressing the velocity potentials 
in terms of the perturbation amplitudes. To find these, consider the kinematic boundary 
condition which states that the normal components of each fiuid's velocity at the interface 
equals the velocity of the interface itself [^. Using the gradient (^ we write the kinematic 
boundary condition for fiow in a sphere as 

1 dp d4)j 



m 

'dt 



sin^ (p/a) dip dip 



p=TZ 



'd_^ 

dp ^ 



(16) 



Inserting expression TZ = R + Ci.Vi't) and equation (^) for 0j into the kinematic boundary 
condition (|16D, we solve for (f)jn{t) consistently to second order in ( to find 



Hn{t) — —-, — f Cn 



QC 



27raS'|?2| 



(17) 



and 



77, 



27ra5'|n| 



+ E ( S^^l?^?^') + rr ) Cn'Cn-n' + 77^^ E ( S^'^l'^^') 



^2 



2C|r2| 



Cn'Cn— n' 



(18) 



The overdot denotes total time derivative. The sign function sgn{nn') = 1 if (nn') > and 
sgn{nn') = —1 if {nn') < 0. 

We can use relations (|T^ and (|T8|) to replace the velocity potentials (pj in the equation of 
motion (|1^) with the perturbation ( and its time derivative (. Keeping only quadratic terms 
in the perturbation amplitude, and equating Fourier modes n on each side of equation (^TJ), 
leads to the differential equation for perturbation amplitudes (n- For n ^ 0, 



n',n') Cn'Cn-n' + G(n,n') Cn'Cn-n' 



(19) 



where 



X{n) 



Q 



27ra252 



{A\n\-C) 



a 



\n\(n^ - 1] 



(20) 



is the linear growth rate, and 
F{n, n') 



\n\ 


\ QAC 


\^ ( '^1 


aC 


aS 


[27ra2^2 


- — sgn[nn ) 





n 



1 - ^(3n' + n) 



G{n, n') = — {A\n\[l - sgn{nn')] - C} 
aS 



+ 



Q 



(21) 
(22) 



are the second-order mode coupling terms. Equation (^) is the mode coupling equation of 
the Saffman- Taylor problem in a spherical Hele-Shaw cell. It gives us the time evolution 

10 



of the perturbation amplitudes (n accurate to second order. In the following sections we 
study equation ([T9|) in more detail, and investigate the role played by cell geometry in the 
interface dynamics. 



III. DISCUSSION 



We use the mode coupling equation ( jT9| ) to investigate the linear instability of individual 
modes and the coupling of a small number of modes. The most noteworthy effect of curvature 
we identify concerns its influence on finger tip-splitting. Tip-splitting is related to the 
influence of a fundamental mode on the growth of its harmonic We take n as the 
fundamental and 2n as the harmonic. To observe interfacial instability of the fundamental 
mode n, we must have A(n) > 0. This occurs if the destabilizing contribution QA in 
Eq. ([20| ) is positive and sufficiently large compared with the stabilizing surface tension term 
proportional to a. To observe growth of the harmonic mode 2n, we presume that QA is 
sufficiently large that X{2n) is non- negative. 

To simplify our discussion it is convenient to rewrite the net perturbation in terms 
of cosine and sine modes 

oo 

C(^, = Co + E Kit) cosM) + bnit) sinine)] , (23) 

n=l 

where ctn = C™ + C-n ^ind bn = i (Cn ~ C-n) are real-valued. Without loss of generality we 
may choose the phase of the fundamental mode so that a„ > and 6„ = 0. We replace the 
time derivative terms a„ and 6„ by X{n) an and X{n) bn, respectively, for consistent second 
order expressions. Under these circumstances the equations of motion become 

dsn = A(2n) a2n + \ T{2n, n) al (24) 



b2n = X{2n) b2n, 
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(25) 



where the tip-sphtting function is defined as 

T(2n, n) = [F{2n, n) + \{n) G{2n, n)] . (26) 

Note that the sign of T{2n,n) dictates if finger tip-sphtting is favored or not by the 
dynamics. If T{2n,n) < 0, at second order the result is a driving term of order forcing 
growth of a2n < 0. With this particular phase of the harmonic forced by the dynamics, 
the n outwards-pointing fingers of the fundamental mode n tend to split. In this case the 
driving term in equation of motion (^) spontaneously generates the harmonic mode. In 
contrast, if T{2n,n) > growth of a2„ > would be favored, leading to outwards-pointing 
finger tip-sharpening. Note that mode h2ni whose growth is uninfiuenced by a„, skews the 
fingers of mode n. In the presence of < 0, the role of is to favor one of the two split 
fingers over the other. 



A. The flat-cell limit 

We begin our discussion by analyzing the fiat-cell limit of the mode-coupling expres- 
sion (|T^). We hold fixed the unperturbed interface velocity v = Q/C. Three distinct fiat 
space limits can be taken: (i) {north pole) let a ^ oo holding R finite so that (7 — 1 and 
aS R; (ii) {south pole) let a — oo holding R' = tto — R finite; so that C ^ —1 and 
aS R'; (iii) {tropical) let R/a be constant and hold fixed k = 2Txn/ C These three limits 
correspond to three physically distinct fiat space fiow problems. Our goal in this section is 
to verify that the linear and nonlinear terms in the equation of motion (^) reduce to their 
expected forms in the fiat space limit. The actual evolution of interfaces according to the 
equations of motion is then discussed in greater detail in sections [III Ij| and p.11 C| . 



At the north pole (case (i), C 1), we recover the mode-coupling equations 0] of fiat, 
radial divergent fiow [^, related to outward radial motion in which fiuid 1 pushes fiuid 2: 
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F{n,n') 



\n\ j QA rl 




(27) 



Provided the viscosity contrast A > 0, so that the less viscous fluid pushes the more vis- 
cous fluid, the interface is hnearly unstable and exhibits flnger growth. As the interface 
perturbation grows, the nonlinear mode coupling broadens and splits the outward-pointing 
fingers and sharpens the inward-pointing fingers. This occurs because T[2n, n) < when 
X{2n) > 0. 

In contrast, at the south pole (case (ii), C —>■ —1), we obtain fiat, radial convergent 
fiow ||2^ equivalent to the inward radial motion corresponding to withdrawal of fiuid 2 
surrounded by fiuid 1. Provided that QA > 0, so that the less viscous fiuid displaces the 
more viscous fiuid, the interface remains linearly unstable. However, now the outwards- 
pointing fingers sharpen because for negative C we find T{2n,n) > when X{2n) > 0. 
The asymmetry between north and south pole behaviors occurs primarily in the nonlinear 
term and comes from the terms proportional to C. In contrast, the linear growth rate 
X{n) is nearly symmetric between the north and south pole limits. This can be understood 
because interchanging the north and south poles is equivalent to reversing the sign of Q 
(the direction of fiow) and the sign of A (interchanging the fiuid viscosities) while holding 
the surface tension a unchanged. However, the term in X{n) proportional to C breaks this 
symmetry slightly. 

A special example of the tropical case (iii) is the equatorial limit, for which C = 
and S = 1. This case should be compared with the cylindrical fiow geometry, taking the 
cylinder tangent to the sphere at the equator. The fiat space limit of this problem indeed 
reduces to the problem of rectangular fiow in fiat space with periodic boundary conditions. 
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In particular, 

(28) 

F{n,n') 
G{n,n') -^A\k\ [1 - sgn{kk')] . 

As expected, there is no tendency for tip-splitting on the cylinder. This alludes to our sug- 
gestion that tip-splitting is controlled by Gaussian curvature K rather than mean curvature 
H because variation of the radius of curvature of the cylinder alters the mean curvature 
while the Gaussian curvature remains zero. Tip-splitting is absent for any value of the mean 
curvature. 

Before the limit is reached, there are small differences between the spherical and cylin- 
drical cases. For example, the surface tension contribution to A(n) vanishes for A; = on 
a cylinder corresponding to translation invariance of the interface length. On a sphere the 
corresponding displacement moves the interface from the equator to the tropics, shortening 
the interface and lowering the surface energy. However, on a sphere the modes n = ±1 
correspond to a global off-center shift of the circular interface preserving circular shape and 
perimeter. Thus the surface tension term vanishes for n = ±1 on the sphere while this mode 
increases the perimeter and raises the energy on a cylinder. 



Av — ak 



B. Linear growth 



Consider the purely linear contribution, which appears as the first term on the right hand 
side of equation (|T^). Since R varies with time, the linear growth rate X{n) is time dependent 
as well. This implies that the actual relaxation or growth of mode n is not proportional to 
the factor exp[X{n)t], but rather 



Cn(t) =Cn(0)exp 



X{n) dt' 



(29) 
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If /q \{n) dt' > the disturbance grows, indicating instability. Two relevant facts can be 
extracted from the linear growth rate: (i) the existence of a series of critical radii Rc{n) 
(defined by setting X{n) = 0) at which the interface becomes unstable for a given mode n; 
(ii) the presence of a fastest growing mode n*, given by the closest integer to the maximum of 
equation (pOD with respect to n (defined by setting d\{n)/dn = 0). In view of equation (p9D 
n* is not simply related to the number of fingers present, even in the early stages of pattern 
formation. Furthermore, in the nonlinear regime the subsequent tip-splitting process and 
mode competition result in a final number of fingers which can differ from the number 
present in the linear regime. 

Cardoso and Woods p4| analyze flat space radial flow assuming the presence of a constant 



low level of noise during the whole evolution of the interface (their "model B", also see 
Ref. M). The sources of noise may come from inhomogeneities on the surface of the Hele- 



Shaw cell, irregularities in the gap thickness b, or from thermal or pressure fluctuations [25 



The predictions of this model are in qualitative agreement with experimental observations 



within the linear regime p4| and the nonlinear regime 0. 

Suppose that we begin with an initially circular interface that is steadily expanding. 
During the interface expansion each mode n is perturbed with a constant (in time) random 
complex amplitude Cn(0)- This noise amplitude contains an n dependent random phase 
but its magnitude |Cn(0)| is independent of n by assumption. As the interface continues to 
expand, it progressively reaches critical radii Rc{n) for n = 2 , 3 , ... . Once a particular 
Rc{n) is reached, the perturbation amplitude (n starts to vary with time. Within this model, 
the first order (linear) solution of equation ([T9|) can be written as 

Cn(0) if R<Rc{n) 

[ Cn(0) {(itj?) [i] exp [C{A\n\ - C) (Ti - 1)]| if i? > i?,(n) 
where the functions Ti = ta.ia[Rc{n)/a]/ taii{R/a), T2 = tan [i?c(n)/2a]/ tan (i?/2a), and 
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Cc = COS [Rc{n)/a]. 

To see the overall effect of Eq. (pOD, we plot the time evolution of the interface using 
the experimental parameters given in Paterson's classical experiment ||^. Paterson observed 
the rapid growth of fingers, as air {rji ^ 0) was blown at a relatively high injection rate, 
Q = 9.3 cm^/s, into glycerine (772 ~ 5.21 g/(cm s)) in a radial source flow Hele-Shaw cell. 
The thickness of the cell b = 0.15 cm and the surface tension a = 63 dyne/cm. We take 
into account modes n ranging from n = 2 up to 20. We evolve from initial radius Rq = 0.05 
cm. The noise amplitude |Cn(0)| = -Rq/IOOO. Figure 2 depicts the evolution of the interface, 
for a random choice of phases, up to time t = 20 s. We set the radius of the sphere a = 10 
cm. We encourage the reader to compare the resulting interface with equivalent figures in 
fiat space. Ref. [0] contains a figure in which growth conditions are similar. In particular, 
the phases of the initial perturbations are identical, although in the present case we choose 
smaller initial amplitudes. It is also of interest to compare with fiat space experimental 
patterns found in the literature [3,^ p3| , p4tp6 . 



C. Nonlinear behavior and tip-splitting 

To visualize the consequences of the second order term in the equations of motion, we 
solve equation (|19D to second order accuracy. If we substitute the linear solution given in 
equation ( pOD into the second-order terms on the right hand side of equation (|19|), we obtain 
the differential equation 

Cn = A(n) C^ + W{n,t), (31) 

where 

W{n, t)=Y: [F{n, n') C^n' + G(n, n') C^"C"-„'] (32) 
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acts as a driving force in the linearized equation of motion (|3T|). Equation (|31|) is a standard 



first order linear differential equation [27| with the solution 



Cn(0) if/?<i?eH 
Cn(t) = <! (33) 



[n 



Here tcl?^) is the time required for the unperturbed growth to reach radius Rc{n) and can be 
calculated from equation (|^) . This solution describes the weakly nonlinear evolution, where 
the dominant modes just become coupled by nonlinear effects. 

We use the second order solution (p3D to investigate the nonlinear coupling among various 
modes n. In figure 3, we plot the interface for a certain time (t = 20 s), considering the 
same random choice of initial phases as was employed in figure 2, and coupling all modes 
with 2 < n < 20. The nonlinear evolution leads to wider fingers and their tips become more 
blunt. These fingers spread and some of them start to bifurcate, by splitting at the tip. As 
was shown in Eq. (0), tip splitting is caused by fingers of mode n driving the growth of 
their own harmonic 2n. 

Now we turn our attention to the investigation of how curvature influences finger tip- 
splitting for Hele-Shaw flow on a sphere. It has been shown in reference that, for fiat, 
divergent radial flow under the condition that X{2n) > 0, we have T{2n, n) < 0. Thus, when 
the harmonic is able to grow, finger tip-splitting is favored. In section [111 A| we showed that 



in the limit of small cell curvature K we recover the fiat space limit where we know T{2n, n) 
is negative. To carry out our current analysis, we consider the case in which \{2n) = and 
see how the curvature influences T{2n,n). 

Since spherical Hele-Shaw flow involves many independent parameters we change only 
one relevant quantity at a time to see what each one does. We hold the unperturbed interface 
velocity v and the unperturbed interface contour length C fixed, while varying the curvature 
K. This isolates the influence of spatial curvature K from the effect of variations in v and 
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L. Fixed L means that fixed mode number n corresponds to a fixed wavelength. 

We adopt an instantaneous approach: we look at the linear growth rate and mode 
coupling at an instant in time, ignoring the past history of how a given interface arose from 
some initial condition followed by growth. Such an instantaneous approach, at which C 
and V have a particular value, enables us to compare the behavior of interfaces evolving in 
distinct background curvatures, but under dynamically equivalent circumstances. Moreover, 
at the instant when the interface circumference C in curved space equals the circumference 
in flat space, if the two velocities also match, the identity Q — vC shows that the Q value in 
flat space equals the Q value in curved space. Therefore, it is advantageous to look at the 
instantaneous tendency towards tip-splitting for interfaces of identical unperturbed v and 
C. 

Consider a particular v and C combination at the onset of growth of mode 2n (using the 
condition A(2n) = 0) in the limit of flat space, where it is known that T(2n,n) < 0. To 
illustrate how tip-splitting varies with curvature for flow on a sphere, we plot in flgure 4 the 
tip-splitting function T(2n, n) as a function of K. The sohd (dashed) curves refer to the 
behavior in the northern (southern) hemisphere. By inspecting flgure 4 we notice that, in the 
northern hemisphere there is a suppression of tip-splitting for increasingly larger curvature 
because the magnitude of T(2n, n) is maximum for the flat, divergent radial case {K = 0) 
and decreases as curvature is increased. For each n there is some value of K for which the 
circle of circumference C hits the equator. This is precisely the point at which northern and 
southern hemisphere branches of the curves meet in flgure 4. 

For the southern hemisphere T{2n,n) is always positive. Consequently, we should not 
expect flnger tip-splitting of outward-pointing flngers in this hemisphere. Actually, flnger tip- 
splitting is replaced by tip-narrowing, along with a splitting of the inward-pointing flngers. 
Finger tip-narrowing is regulated by the curvature K. 
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Another noteworthy point about figure 4 is the evident symmetry-breaking in T{2n, n) 
between northern and southern hemispheres. The justification for this asymmetry is similar 
to the one that explains why the interface is always in-out asymmetric in fiat radial fiow: if 
we are located at either pole we can always distinguish regions that are inside and outside 



the interface. That is why convergent and divergent fiat, radial fiows are not equivalent [23 
From this point of view, the asymmetry observed in figure 4 should be somehow expected. 

We conclude this section presenting the lowest order curvature expansion for T{2n, n) 
and studying how curvature K infiuences tip-splitting for small curvature values. Using the 
instantaneous approach described above, series expansion of equation ( ^6]) for the northern 
hemisphere yields 



I {2n,n) = ; -r 

^ ' ^ 32iT^a n{2n + 



Svr^a n2(2n + 1)2 



+ K 



+ 0{K^). (34) 



The perturbative expansion (p^ ) for small K explicitly shows a linear correction to the fiat 
space limit {K = 0) that reduces the magnitude of T{2n,n), inhibiting finger tip-splitting. 
A similar kind of expansion can be done for the southern hemisphere, resulting in enhanced 
tip-sharpening. Those results are in agreement with the small curvature behavior depicted 
in figure 4. 

Equation ( P^ allows us to make predictions about tip-splitting behavior for Hele-Shaw 
fiow in (weakly) negatively curved backgrounds. This type of fiow could happen, for instance, 
between two saddle-like surfaces. Actually, fiow in porous media seems to be somehow linked 
to fiow on negatively curved surfaces. Porous materials define multiply-connected surfaces, 
presenting negative average Gausssian curvature 

JKdS An (1 - g) 

where dS denotes an infinitesimal area element and g is the so-called genus [|r^, which 
denotes the number of holes present in a given surface. Expression (|35|) relates the integral 

19 



of the Gaussian curvature K of a given surface to its topological properties, by virtue of 



the Gauss-Bonnet theorem [19|. From (|35D we verify that K becomes progressively negative 
when the number of holes is increased: that is why ^ > for a sphere {g = 0), K = 
for a torus {g = 1) and ^ < for a (7-torus {g > 2). In this sense, a medium which is 
rich in pores (holes) would present negative curvature features. However, flow on surfaces of 
constant negative Gaussian curvature is complex and to treat the problem rigorously would 
require an interesting generalization of our formalism beyond the scope of the present paper. 
Here we simply point out what would be the behavior for flow on surfaces that are slightly 



(negatively) curved. It is easy to see, by performing the substitution K — > —K in (0), that 
negative curvature should enhance finger splitting in comparison to the flat case. This is 
in striking contrast to the problem of flow in positively curved surfaces, such as the sphere, 
where curvature leads to suppression of tip-splitting. This sets an important distinction 
between Hele-Shaw flows on surfaces with negative curvature and those on surfaces with 
positive K. Our results confirm explicitly the predictions made in reference about the 
role of cell curvature on finger splitting. 

The fact that the correction is linear in K supports our suggestion that Gaussian curva- 
ture is more relevant than mean curvature. The same expansion Eq. (|3^) could be written 
in terms of the square of mean curvature H. However, then we would predict supression of 
tip-splitting for both positive {H > 0) and negatively {H < 0) curved surfaces. We believe 
tip-splitting is enhanced for negative Gaussian curvature because the metric creates ample 
space for the spht fingers to penetrate without mutual competition. 

IV. CONCLUDING REMARKS 

In this paper we generalized the traditional Saffman- Taylor problem by studying vis- 
cous flow on curved surfaces. Our main purpose was to investigate the influence of spatial 
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curvature on viscous fingering pattern formation, wlien fluid flow takes place on a sphere. 
By deriving the equation of motion for the interface pertubation amphtudes, using a mode- 
coupling approach, a study of both linear and weakly nonlinear stages of evolution could be 
carried out. 

We have shown that cell curvature can be used as a control parameter to discipline 
splitting of the viscous fingers. The fluid-fluid interface can be more stable or unstable, 
with respect to tip-splitting, depending on the curvature of the surface to which the flow 
is confined. We also detected an asymmetry on tip-splitting behavior depending where 
the interface evolves: while tip-splitting may be still present on the northern hemisphere, 
it is completely replaced by finger tip-sharpening on the southern hemisphere. We also 
found evidence that Hele-Shaw flows on negatively curved surfaces would present enhanced 
tendency to tip-splitting, so highly branched patterns may be expected. In summary, we 
have explicitly verified that interfacial behavior is coupled to the geometry of the Hele-Shaw 
cell, so that curvature has important consequences for flow dynamics. 
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APPENDIX A: DARCY'S LAW ON A SPHERE 



This appendix derives Darcy's law for viscous flow between concentric spheres. The 



derivation is based upon results presented in Bird, et al. |28]. We begin with a coordinate- 



free representation of the continuity equation for an incompressible fluid 

V ■ u = (Al) 
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and the Navier-Stokes equation 

-Vp + ?7V^u (A2) 



P 



•9u , 
g^ + (u.V)u 



where u denotes the three-dimensional fluid velocity and we neglect the acceleration due 
to gravity. Neglecting the inertial terms on the left-hand side of Eq. (^2), transforming to 



spherical coordinates (r, 6, (p), and specializing to the case of polar flow {ur = u^^ = 0) p8| , |29 
we rewrite the Navier-Stokes equation 

Zl^ = r]\^ ue- ^ . 2n ]- (A3) 



r do V r2sin^6'. 
The noteworthy term in Eq. is the second term multiplying viscosity which enters as 

a result of the curvilinear coordinate system. 

To determine U0{r,9) note that the continuity equation ([AID adapted for polar flow 

1 d{ug sin 9) 



r sin 6 89 

demands a solution of the form 



(A4) 



= (A5) 

Insert this form into the Navier-Stokes equation ( |A3| ) and multiply by rsin6' to separate 
variables: 

d9 r dr^ dr^ ^ ^ 

Because the left-hand side involves only angle 9 and the right-hand side is radial, involving 
only r, each side must be constant sharing a common value B. The solution of the radial 
equation subject to no-slip boundary conditions u{a) = u{a + 6) = is 

uir) = ^ir-a)ia + b-r) ^^^^ 
b r 

where B = —2riauQ/b'^. 
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Darcy's law for Hele-Shaw cells is obtained by averaging the three-dimensional velocity u 
over the gap width b. We choose to average in the radial direction r and obtain vq = u/ sm9 
with 



6 [a 6 V 2a 



(A8) 



where 



2a 



log 1 



a , 



(A9) 



Equation ( [X9| ) could also be written, in a more compact form, as a hypergeometric function 
J^{b/a) =F(l,2;4;-6/a). 

Equation (^) generalizes the usual flat-cell velocity average. Darcy's law becomes 



127] \ade 



(AlO) 



The effect of curvature can thus be incorporated entirely into a reduced gap width b or an 
enhanced viscosity r/. Equation ( |A10| ) recovers the usual Darcy's law for flow in flat Hele- 
Shaw cells with corrections of higher order in b/a. In this work we will be interested in the 
case b <^ ( <ti p, with p a. 
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FIGURE CAPTIONS 



FIG. 1: Schematic configuration of tlie flow in a splierical Hele-Sfiaw cell. 

FIG. 2: Time evolution of the fluids according to Eq. (|30|), including modes 2 < n < 20. 
The initial perturbation amplitudes |Cn(0)| = -Ro/1000 and Ro = 0.05 cm. Other parameters 
are given in the text. We show the fluid- fluid interface at t = 20 s. 

FIG. 3: Nonlinear evolution according to Eq. (pSj). All physical parameters and initial 
conditions are the same as those used in Fig. 2. 

FIG. 4: Variation of T{2n,n) as a function of the spherical Hele-Shaw cell curvature K, for 
modes n = 8, 10, 12. The solid (dashed) curves describe behavior in the Northern (Southern) 
hemisphere. The units of T{2n,n) and K are (cm s)~^ and cm~^, respectively. 
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